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ECOMB is established as a hadronic multiparticle production generator by soft 
interaction. It incorporates the eikonal formalism, parton model, color mutation, 
branching, resonance production and decay. A partonic cluster, being color-neutral 
initially, splits into smaller color-neutral clusters successively due to the color mu- 
tation of the quarks. The process stops at hadronic resonance, qq pair, formation. 
The model contains self-similar dynamics and exhibits scaling behavior in the fac- 
torial moments, e.g. the intermittency. 



1 Introduction 

The study of multiparticle production in low-p^ processes has been pursued for 
over twenty years. Since they involve soft interactions, they cannot be treated 
in perturbative QCD. In the absence of any reliable theoretical approach to 
the problem, many models have been proposed, .most of which are represented 
in the review volumes published ten years ago El. Nearly all of those models 
have since been shown to be inadequate in light of the data on fluctuations and 
intermittency □. Indeed, there are very few models that have the appropriate 
dynamical content capable of reproducing the scaling behaviors observed in 
the experiments. In this paper we propose a model that incorporates some 
aspects of non-perturbative QCD and is capable of generating the features of 
the intermittency data, which are shown in the last figure of this paper. To our 
knowledge those data have not been fitted by any model that contains some 
features of the color dynamics. To reproduce those data has become the main 



motivator for this work. 

Our approach u embraces many time-honored properties of hadronic col- 
lisions. Since hadrons are extended objects, the eikonal formalism is at the 
foundation of our model. Thus it is not difficult for our model to possess 
the virtues of geometrical scaling, approximate KNO scaling and unitarity. In 
order to build into the model features of chromodynamics, it is necessary to 
introduce quarks and gluons into the eikonal formalism, so the parton model is 
an essential gateway into the microscopic domain of color interactions. Once 
we enter that domain, we embark on an unconventional journey of studying 
color mutation of the constituents as a dynamical process by which the col- 
ors of the quarks evolve through the emission and absorption of gluons. The 
smallness of as is never assumed, so the evolution is not perturbative. With 
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all partons taken into consideration globally, we follow the evolution of the 
configuration in the color space. When the configuration exhibits color neutral 
subclusters, we allow branching to take place with a possible contraction of 
the cluster size in accordance to a reasonable rule consistent with confinement 
dynamics. Successive branching leads to smaller and smaller clusters, until 
they are finally identified as particles and resonances. The decay of resonances 
are also taken into account before the final state of an event is determined. 
Evidently, the model contains many features of soft interaction that are famil- 
iar and desirable at a qualitative level. Here we put them on a quantitative 
basis. 

We shall call this model ECOMbI, which stands for cikonal color mutation 
branching. Since the partons play a fundamental role in this model, it is 
significantly closer to QCD than any eikonal model on soft hadronic collisions 
has ever been. 

There are a few parameters in the model. They are adjusted to fit a 
large body of experimental data on low-pr processes for y/s ~ 100 GeV. They 
include a e i, UineU ( n )i C q , dn/dy, P ra , and F q for all s in the range 10 < \fs < 
100 GeV, and all rapidity intervals. F q are the normalized factorial moments, 
whose power-law dependence on the rapidity bin-size 5 has been referred to as 
the intermittency behavior □. Except for F q , all the other pieces of data are 
global in nature; i.e., examined in or averaged over all rapidity space. They 
can be fitted by many models. F q in small rapidity intervals exhibit local 
fluctuations, which are what invalidate most of those models. 

2 Geometrical Effect on pp Soft Collision 

To focus on processes in which hard subprocesses are unimportant, we con- 
fine our attention to the energy range 10 < \fs < 100 GeV, which covers the 
CERN ISR energies. From parton model point of view, this soft process with 
a relative smaller momentum transfer undergoes a relative longer time evolu- 
tion of partons than hard subprocess. High energy hadronic collision provides 
a setting up of the initial condition for the partonic evolution. A proper de- 
scription of the collision is the basis of any soft interaction model. Two of the 
well known experimental results, geometrical scaling and approximate Koba- 
Neilsen-Olesen (KNO) scaling, provide immediate constrains on the collision 
model building. The eikonal formalism of hadronic collisions gives us an ideal 
geometrical framework to fit all the constrains. Here is the highlight of it. 

In terms of the eikonal function £l(b), the geometrical scaling, i.e., u e i/atot 
is roughly constant, can be guaranteed, if fl depends only on the scaled impact 



2 



parameter R — b/bo(s). The reduced inelasticity function 



g(R) = 1 - e 



(1) 



satisfies the normalization condition J °° dR 2 g(R) — 1 by setting Omei = 
7t6q(s). The function g(R) describes the probability of having an inelastic col- 
lision at R. For pp collisions the well-determined form for the eikonal function 



which has been used to give a good description of da/dtu. 

The /ith term of the inelasticity function expansion in a power series, 



may be regarded as thjeatth-order rescattering contribution, i-and can be related 
to the /j,-cut-PomeronH0. The average /i is estimated as 1.60, according to (|J). 

To go further in multiparticle production, it is necessary to model the 
dynamics of soft interaction. This paper is devoted to develop such a model 
by incorporating the parton model into the formalism as far as we can. 

To do so, we adopt the eikonal description to refer to parton number; let 
B£ be the probability of having v partons in the /u-cut Pomeron when the two 
incident hadrons are separated by a scaled impact parameter R. If we use 
£v— >n{" • •} to denote the evolution process that takes the v partons to the n 
hadrons, we may express symbolically the multiplicity distribution of hadron 
as 



The form invoked in the bracket essentially provides a parton number distri- 
bution which involved in the /i-cut Pomeron exchange soft interaction and the 
description of £ v ^ n is the main task of this paper. 

A way of thinking about B£ is to consider the /i = 1 case, for which we 
have an one-Pomeron exchange diagram for the elastic scattering amplitude 
with the Pomeron being cut to expose the internal quark, antiquark and gluon 
lines on mass shell. We use B£ to represent the probability having v such 
partons as fi Pomerons in the elastic amplitude are cut. 

Since there exists no rigorous derivation of Bp, we shall assume that it is 
Poisson distributed around some mean number v of partons. It is reasonable, 




(2) 




(3) 
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since it is known that a cut ladder corresponds to a multiperipheral diagram for 
a ^-particle production amplitude, whose rapidity distribution is uniform, and 
multiplicity distribution Poissonian. The mean number can depend on both \x 
and s. Being guided by so-called logs physics, we adopt a parameterization for 

P(/i,*) = *(»KW , ( 5 ) 

where the parameters depend on the energy by v(s) = isq + i^lns and a(s) = 
ao + ailns + a2ln 2 s. 

Once the number of parton, v, is determined, one can move on the issue 
of evolution of the v parton cluster. Before going further, I like to point out 
an important departure from the usual application of the parton model: We 
combine the property of hadrons being spatially extended objects with the 
property of hadrons being made up of smaller constituents. Clearly, for most 
collisions whose impact parameters are nonzero, only portions of the partons 
in the incident hadrons overlap and interact, so the number of hadrons pro- 
duced depends on that overlap. The event-to-event fluctuation of the particle 
multiplicity therefore depends strongly on the impact-parameter fluctuation, 
and the event averaged parton distributions, as determined from the structure 
functions, are of no use in that respect. 

As one of the consequences of the carefully consideration of event-to-event 
fluctuation, not only the average charged multiplicity of hadron (n c h), but also 
the KNO scaling are reproduced. They are satisfactorily consistent with the 
datacl in terms of the standard moments C q (s) = (n q ) / (n) q over an energy 
range of 10 < \fs < 70 GeV as shown in Fig 1. 

3 Color Mutation 

We now consider the color dynamics of soft interaction. The initial state is 
that there are v partons distributed in some fashion in a linear array in rapidity 
space. This is a consequence of the parton model in that the partons are in 
the incident hadrons to begin with, and the collision rearranges those partons 
and sets off the evolution process that takes those partons to the final state, 
where the produced hadrons are decoupled. The evolution process involves 
quarks and antiquarks emitting and absorbing gluons. Since the process is 
not perturbative, there is no analytical method to track the time development 
of the process. Thus we shall use Monte Carlo simulation to generate the 
configuration at each time step. 

There are two spaces in which we must track the motion of the color 
charges. One is the two-dimensional color space; the other is the one-dimensional 
rapidity space. The latter can be extended to include the azimuthal angle <j> 
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and the transverse momentum \pr\, but in our first attempt here we integrate 
over those variables and examine the simpler problem of a ID system. As for 
the color space, it is 2D, and spanned by the (h,Y) axes, as in SU(3) flavor. 

The configuration in the two spaces are coordinated in the following sense. 
Starting from the extreme left end of the rapidity (77) space, for which the 
system that contains no partons is by definition color neutral and therefore is 
represented by a point at the origin of the color space, we move in the positive 
direction in 77 space until we cross a color charge. At the point to the right 
of that color charge, the corresponding point in the color space jumps from 
the origin to a position that represents the color of the charge that has been 
crossed. Let us denote that jump by a vector. Then each time we cross a color 
charge in the 77 space, there is a corresponding vector added to the previous 
point in the color space. The succession of additions of those vectors, each one 
starting from the tip of the previous one, forms a path. Since the whole parton 
system is color neutral, by the time we have moved to the extreme right in the 
77 space, the path in the color space returns to the origin, thus forming a closed 
loop. Such a loop may be self-intersecting at various points. We call such a 
closed path a configuration of the system in color space. 

In order to initiate the time evolution, we need distribute the v partons 
into both of those two spaces to form such an initial configuration. 

(1) Initial Distribution in Rapidity Space 

By rapidity we mean space-time rapidity 77, where 77 = \f-'njzr, because we 
are dealing with a variation of spatial separations between partons due to the 
color force during the time evolution. 

The usual parton distribution in momentum fraction (at the lowest acces- 
sible Q 2 ) implies roughly a flat 77 distribution with rapid damping toward zero 
at large 77. It is modeled as 



p(v)=P {1) (v)+P {2) (r ] ), 



(6) 



P W (v) = Po( 



77 + 770 
2770 



% + 770) 6>(?7o - 77) + % - 770) % c - 77) + 



Tjmax 7] 



0(77 - Tj c )e(i] max - 



v)), 



(7) 



Tjmax 1]c 



P (2) (v)=Po( 



-V + Vo e 
2??o 

ffmax H~ T) 



8(V + Vo) %o -v) + 0(V + Vc) 9(-Vo - v) + 



0{rj + r] max ) 9{-rj c - 7/)) . 



(8) 



Vmax 
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in which, (77) and p^ (77) indicate the forward and backward semi-spheres 
respectively. Only small portion of partons overlap in a limited central region 
(—770,770) after collision follows from the known fact that the correlation be- 
tween produced particles is short-ranged, because the soft interaction range is 
short in the conventional wisdom. 

In the color space spanned by 73 and Y , a quark is represented by a vector, 
which has the coordinates of one of the triplets: (1/2, 1/3), (—1/2, 1/3), and 
(0, —2/3). An antiquark is represented by a vector directed opposite to one of 
the above. As a distribution of partons in the 77 space is generated, say, p^ (77) 
between — 77 max < V < Vo, we assign to each parton a color vector consistent 
with the requirement that the quarks and antiquarks come in pairs, but their 
orderings in the 77 and the color spaces are totally random. The partons in 
p^{rf) are distributed similarly, but completely independently, between —770 < 
t? < 77 max . Due to the merge in the central region (—770,770), v partons are no 
longer divided into two separate color neutral groups. Usually, they mix up 
into a whole neutral system according to the rule to count the color charge as 
described in the beginning of this section. 



(3) Color Interaction 

We now consider the evolution of a configuration due to QCD dynamics. A 
nonperturbative treatment of v simultaneously interacting color charges is, of 
course, too difficult to contemplate here. We reduce the problem by considering 
pairwise near-neighbor interaction via the exchange of a gluon in any of the 
s-, t-, or 7i-channel, whichever is applicable. Starting from the extreme left in 
the 77 space, we regard the ordered chain of v partons as having v — 1 links 
(with varying link lengths). Pairwise near- neighbor interaction means that we 
consider the v — 1 links in 77 space one at a time, according to a rule to be 
specified below. After the interactions at all v — 1 links are considered, the 
evolution of the whole configuration is regarded as having taken one time step, 
and the process is then repeated. 

For every global configuration a. of the v partons, there is an associated 
energy E a . E a is determined by 



i=l 



i = l 



(9) 



where Cj is the color vector for jth parton, Ci is the color charge the i-link 
seeing on its left side and sums obviously depends on the particular path 
in color space. Whenever a local pair of partons interact, the outcome may 
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or may not affect the global configuration. Our statistical approach to the 
determination of the global configuration a, consistent with favoring the lowest 
energy state, requires that the probability for configuration a to occur is 



where (3 is a free parameter. Using the Metropolis algorithm, we use P a to 
determine the outcome of a local interaction at every link. 

(4) Spatial Fluctuations 

In the meanwhile the system undergoing color mutation due to the exchanging 
of gluon, the spatial position of quarks also fluctuates due to the same reason. 
We track the change in the length of the ith link in the r\ space, di = r/i+i — r/i, 
to quantify the spatial evolution of ith link. Whether the link length contracts 
or expands depends on the attractiveness or repulsiveness of the net color 
forces that act on the two ends of the link. To determine the nature of that 
force is beyond the scope of this treatment. We shall model the change in 
link length by a stochastic approach, consistent with how we have handled the 
color mutation part of the dynamics of the complex system. We allow di to 
change by m;, i.e., di — ► di + mi, where the probability for rm is specified by 
a distribution 



V i {m)= 11 6{m + d i )e{ 1 2-m) + [l- ll { l2 +d i )]5 mi - di , (11) 



which satisfies J2m=-d - ^i( m ) = 1- This distribution allows m to be positive 
(expansion) uniformly up to 72, and negative (contraction) down to — rf,. 71 
and 72 are two free parameters. 

4 Branching and Hadronization 

A fission of the color-singlet cluster occurs, when it contains two color-singlet 
subclusters, since no confining force exists between them. Thereafter, the 
color mutation process is applied to the two subsystems separately and inde- 
pendently. This is repeated again and again until all subclusters consist of 
only qq pairs. At every stage of the evolution process, a cluster shrinks due to 
overall spatial contraction. We always keep the center of the cluster invariant 
to conserve momentum. When a branching occurs, the two daughter clusters 
will have their own respective centers, which will remain invariant during con- 
traction, until they themselves branch. At the end when a hadron is formed 



c 





(10) 
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from a qq pair, the hadron momentum rapidity y will be identified with a value 
of 77 taken randomly between the 77 values of the quark and antiquark. The 
reason for doing this is that without knowing the mass and the transverse mo- 
mentum pt of the hadron. A rough identification of y with 77 is justified for 
free particles at high energy. 

When branching process terminates, the qq pairs are identified as hadrons, 
which may be pions, kaons, or resonances. Those resonances must be allowed 
to decay before the total number and distribution of particles are counted for 
the final state of the event. The probabilities of producing various resonances 
and stable (in strong interaction) particles have been studied experimentally in 
pp collisions intX We use that reference as a generic guide for the proportions of 
all particles produced in any general hadronic collision. We give each hadron 
a transverse momentum according to a "standard" exponential distribution 
of transverse momentum squared. The average transverse momentum is a 
parameter to be fixed actually as 400MeV. If the qq pair forms a resonance, then 
we assume an isotropic decay distribution in the rest frame of the resonance. 
The azimuthal angle is assigned randomly. After boosting back to the cm 
system, the 3-momenta p of the decay particles are then determined. 

5 Results 

In the previous section we have introduced six parameters: r/ c , ?y max and 770 in 
(0) and (g) , (3 in ( |l0| ) , 71 and 72 in ( |ll| ) . They are to be varied to fit the data 
on inclusive distributions and on fluctuations of the exclusive distributions. 

The parameters ?7 C , r7 max and (fcy) = 400MeV are essentially kinematical; 
they set the boundaries of the phase space in which the partons are placed. 
They do not affect the dynamics of color mutation and the spatial fluctuation 
of the phisters. The data used to determine them are the rapidity distribution 
rln/rlrfi 1 ^ anH the transverse momentum distribution. The energy range of the 
data inLL9 is 22 < ^/s < 63 GeV, for which hard scattering is negligible. The 
values of r\ c and ?7 m ax depend on s. In Fig. 2 we show the pseudorapidity dis- 
tributions to compare with ECOMB as an example. The agreement is clearly 
satisfactory. 

The parameters that influence the evolution process of mutation and branch- 
ing are r] , (3, 71 and 72. t]q specifies the overlap region of the initial color- 
neutral clusters, /3 pertains to the probability of color mutation, and 71 and 72 
characterize the fluctuations of the link lengths. The data used tcL determine 
them are on factorial moments F q for q = 1, 2, ... and their variantsa. It is here 
that we can underline the importance of those data on fluctuations, without 
which we have no guidance on how to restrict the detail dynamics of particle 
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production. Putting that in another way, in the absence of a procedure to 
calculate from first principles, any model that fails to fit the fluctuation data 
is missing some aspect of the basic dynamics. To the extent of our awareness, 
very few models on soft interaction have been put to the test of confronting 
those data on the factorial moments for varying bin sizes. 

The parameters that we have varied to give the best fit of the intermittency 
data at y/s = 22 GeV are: v{s) — 9.1, a(s) = 0.63,?y max = 5, r) c — 3.5, 770 = 
1.9, (3 = 0.0015, 71 = 0.077, 72 = 5. In_Fig. 3 we show the intermittency results 
calculated with ECOMB and data ofEJ "which to our knowledge have not been 
reproduced by any model, except ECCOLm. It should be noted that to achieve 
the fits attained is highly nontrivial. If any part of the dynamical process in 
generating the hadrons is altered, one would not be able to obtain the rising 
factorial moments, no matter how many parameters are used. By working 
with the many parts of our model, all of which affect the determination of F q , 
we have gained confidence in regarding the dynamics of color mutation and 
branching as having captured the essential properties of soft interaction. 

6 Conclusion 

It is rather satisfying that we have been able to reproduce the intermittency 
data in Fig. 3. Scaling behavior of that type implies self-similarity in the 
dynamics of particle production. 

We have amalgamated many concepts that form various elements of the 
conventional wisdom about soft interaction. They include: (a) hadrons having 
sizes, (b) eikonalism, (c) parton model, (d) interaction of quarks via gluons, 
(e) statistical properties of a many-body system, (f) spatial contraction and 
expansion of a color system, and (g) resonance production. They are interlaced 
by intricate connections described in this paper. 

The most important omission in this paper is Bose-Einstein correlation, 
without which our model is not completely realistic. That defect must be 
corrected in an improved version. 
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Figure 3: Intermittency data of normalized factorial moments for q = 2 — 5. The data are 
from Ref. [11]. The lines are determined from ECOMB. 
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